function bnet = mk_cancer()
N=3;
dag = zeros(3,3);
dag(1,2)=1;
dag(2,3)=1;
node_sizes= zeros(1,3);
node_sizes(1)=2;
node_sizes(2)=2;
node_sizes(3)=2;
bnet = mk_bnet(dag, node_sizes);
CPT= zeros(2,1);
CPT(1)=0.75;
CPT(2)=0.25;
CPT=reshape(CPT, 1, prod(size(CPT)));
bnet.CPD{1} = tabular_CPD(bnet, 1, 'CPT', CPT);
CPT= zeros(2, 2);
CPT(2,1)=0.12;
CPT(2,2)=0.88;
CPT(1,2)=0.36;
CPT(1,1)=0.64;
CPT=reshape(CPT, 1, prod(size(CPT)));
bnet.CPD{2} = tabular_CPD(bnet, 2, 'CPT', CPT);
CPT= zeros(2, 2);
CPT(2,1)=0.20;
CPT(1,2)=0.76;
CPT(1,1)=0.24;
CPT(2,2)=0.80;
CPT=reshape(CPT, 1, prod(size(CPT)));
bnet.CPD{3} = tabular_CPD(bnet, 3, 'CPT', CPT);
